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ABSTRACT 

Context. Mass-loss occurring in red supergiants (RSGs) is a major contributor to the enrichment of the interstellar medium in dust 
and molecules. The physical mechanism of this mass loss is however relatively poorly known. Betelgeuse is the nearest RSG, and as 
such a prime object for high angular resolution observations of its surface (by interferometry) and close circumstellar environment. 
Aims. The goal of our program is to understand how the material expelled from Betelgeuse is transported from its surface to the 
interstellar medium, and how it evolves chemically in this process. 

Methods. We obtained diffraction-limited images of Betelgeuse and a calibrator (Aldebaran) in six filters in the N band (7.76 to 
12.81 jim) and two filters in the Q band (17.65 and 19.50//m), using the VLT/VISIR instrument. 

Results. Our images show a bright, extended and complex circumstellar envelope at all wavelengths. It is particularly prominent 
longwards of ^9-lOyum, pointing at the presence of 0-rich dust, such as silicates or alumina. A partial circular shell is observed 
between 0.5 and 1.0'' from the star, and could correspond to the inner radius of the dust envelope. Several knots and filamentary 
structures are identified in the nebula. One of the knots, located at a distance of 0.9'' west of the star, is particularly bright and 
compact. 

Conclusions. The circumstellar envelope around Betelgeuse extends at least up to several tens of stellar radii. Its relatively high 
degree of dumpiness indicates an inhomogeneous spatial distribution of the material lost by the star. Its extension corresponds to an 
important intermediate scale, where most of the dust is probably formed, between the hot and compact gaseous envelope observed 
previously in the near infrared and the interstellar medium. 

Key words. Stars: individual: Betelgeuse; Stars: imaging; Stars: supergiants; Stars: circumstellar matter; Stars: mass-loss; Infrared: 
Stars 



1. Introduction 

Several mysteries remain to be solved for evolved supergiant 
stars, among which the structure of their convection, the mecha- 
nism of their mass loss (Levesque 2010 ) and the way dust forms 
in their environment (Verhoelst et al. 2009 ). Josselin et al. ( 2000 ) 
showed that the dust mass-loss rate is not correlated with lu- 
minosity, and that the molecular gas-to-dust ratio shows a very 
large scatter, which is higher than what is observed in asymptotic 
giant branch (AGE) stars. The mass loss mechanism of AGE 
stars, based on pulsations and radiation pressure on dust grains, 
is unlikely to be applicable to red supergiant (RSG) stars because 
they are irregular variables with small amplitudes. The mass loss 
mechanism in RSGs could be related to the turbulent pressure 
generated by convective motions, which, combined with radia- 
tive pressure on molecular lines, may levitate the gas (Josselin & 
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Plezl 2007T) . The red supergiant Betelgeuse (spectral type M1-M2 
la-Ib, Teff = 3650 K, Levesque et al. 2005 ) is known to have a pe- 
culiar convective pattern on its surface characterized by small to 
large convective cells (up to more than one stellar radius in size) 
with supersonic velocities and shocks (Chiavassa et al. 12010b . 
Moreover, the magnetic field discovered on Betelgeuse (Auriere 
et al. 2010 ) may also contribute to the mass loss through Alfven 
waves (see e.g. Hartmann et al. 1984, Pijpers et al. 119891 Cuntz 
et al. 1997, Airapetian et al. 2010). Progress on our understand- 
ing of this intricate combination of physical phenomena requires 
high angular resolution studies of the surface of the star using 
long-baseline interferometry (e.g. Haubois et al. |2009j), and its 
circumstellar environment (hereafter CSE). 

Evidence for mass-loss is well established in late-type stars 
with dust detected tens and hundreds of stellar radii away from 
the stellar surface. However, the way dust is created still re- 
mains a partial mystery, although progress has been made re- 
cently on both the observational and theoretical sides. The 
concept of MOLsphere was first hypothesized spectroscopi- 
cally in the environment of the red supergiants Betelgeuse and 
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[I Cep (Tsuji l20Q0b . Spectro-Interferometric observations con- 
firmed this hypothesis and estabUshed some characteristics of 
the MOLsphere such as its height and composition (Perrin 
et al. 12004, Verhoelst et al. 2006, Perrin et al. I2Q07I see also 
Tsuii 12006 and Schuller et al. 2004) . This scenario is a candi- 
date to explain how dust forms in red supergiants, but does not 
provide rock- solid arguments to explain how the MOLsphere 
can be formed and levitated above the supergiant photosphere. 
One hypothesis is that convection cells detected by Chiavassa 
et al. (.2010) are at the root of mass transfer from the star sur- 
face to its upper atmosphere (Ohnaka et al. 2009 ), and that the 
expelled material condenses as dust within the MOLsphere and 
above. 

In January 2009, Kervella et al. (l2009l hereafter "Paper I") 
obtained adaptive optics (AO) images of Betelgeuse with 
NACO, taking advantage of the cube mode of the CONICA cam- 
era to record separately a large number of short-exposure frames. 
With a selective imaging approach, they obtained diff'raction- 
limited images over the spectral range 1.04 - 2.17yum in 10 
narrow-band filters. These AO images show that the inner CSE 
of Betelgeuse has a complex and irregular structure, with in par- 
ticular a bright "plume" extending in the southwestern quadrant 
up to a radius of at leasy six times the photosphere. These obser- 
vations allowed the probable identification of the CN molecule, 
which is relatively common in the envelopes of evolved post- 
RGB stars (e.g. Cherchneff' 2006 ). We have here a first piece of 
evidence for the mass-loss puzzle of Betelgeuse. A next step re- 
quires detection of the dust itself, as an indicator of the location 
of formation. Considering the temperature at which dust is sup- 
posed to form, this requires thermal infrared imaging. This is 
the object of the new VISIR observations in the N and Q bands 
reported in the present article. The observations and raw data 
processing are presented in Sect. [2] their analysis in Sect. [3] and 
we discuss our results in Sect.lU 



2. Observations 

2.1. Instrumental setup and data acquisition 

For our observations of Betelgeuse (aOri, HD 39801, HR 
2061) and its point spread function (hereafter PSF) calibrator 
Aldebaran (aTau, HD 29139, HR 1457), we used the VISIR in- 
strument (Lagage et al. 2004), installed at the Cassegrain focus 
of the Melipal telescope (UT3) of ESO's Very Large Telescope 
(Paranal, Chile). VISIR is a mid-infrared (hereafter MIR) im- 
ager and slit spectrometer, that can in principle reach a very high 
angular resolution, thanks to the 8 m diameter of the telescope. 
However, under standard conditions at Paranal (median seeing 
of 0.8'' at 0.5yum), the 8 m telescope is not diff'raction limited 
in the MIR (seeing ^ 0.4'' vs. 0.3" diff'raction). A few mov- 
ing speckles and tip-tilt usually degrade the quality of the image 
(see e.g. Tokovinin, Sarazin & Smette 2007 ). To overcome this 
limitation, a specific mode of the instrument, the BURST mode 
(Doucet et al. 2007 ), has been implemented to give the possibil- 
ity to record separately a large number of very short exposures 
(At < 50 m s), in order to freeze the turbulence (see e.g. Law 
et al. 120061 Law et al. 120091 and references therein). To obtain 
diff'raction limited images of the CSE of Betelgeuse, we there- 
fore selected this instrument mode. 

We observed Betelgeuse and Aldebaran on 12 November 
2010, between UT 06:54 and 08:04 (Table [T]). We recorded 
BURST image cubes of several thousand images in eight VISIR 
filters: J7.9 (Ao = 7.76yum, half-band width = 0.55 yum), 
PAHl (^0 = 8.59 ± 0.42 yum), SIV (^o = 10.49 yum, A^ = 



Table 1. Log of the VISIR observations of Betelgeuse and 
Aldebaran. 
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Notes. MJD* is the modified Julian date of the start the exposures on the 
target minus 55 500. The Detector Integration Time (DIT) is the expo- 
sure time for one BURST image, 6 the seeing in the visible {A - 0.5 jjm) 
measured by the observatory DIMM sensor, and AM the airmass of the 
observation. 
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Fig. 1. Transmission of the VISIR filters used for the observa- 
tions of Betelgeuse. The transmission of the J7.9 filter is sup- 
posed to be Gaussian with a central wavelength of 7.76 yum and 
a half-band width of 0.55 yum. The ISO spectrum of Betelgeuse 
is also shown for reference, with an arbitrary intensity scale (thin 
grey curve, see Sect. |3.2| for details). The central wavelength of 
each filter is given in the legend. 



0.16yum), PAH2 (^o = 11.25yum, A^ = 0.59 yum), SIC (^o = 
11.85yum, AX = 2.34 yum), Nell (.lo = 12.81 yum, A^ = 0.21 yum), 
Ql (.lo = 17.65 yum, A^ = 0.83 yum), and Q3 (^o = 19.50yum, 
AA = 0.40 yum). Their transmission curves are shown in 
and are available from the VISIR instrument web site at 
The images were chopped and nodded in a North-South and 
East- West direction by off'setting the M2 mirror and the tele- 
scope (respectively) in order to remove the fluctuating thermal 
background in the post-processing. The chopping and nodding 
amplitudes were both set to 8" on the sky. The chopping fre- 
quency was set to 0.5 Hz in the N band and 1 Hz in the Q band, 
and the nodding period was set to 90 s. 



http ://w w w. eso. org/sci/f acilities/paranal/instruments/visir/ 



2 



p. Kervella et al.: The close circumstellar environment of Betelgeuse from 7.76 to 19.50yL/m 



2.2. RssN data processing 

The data processing procedure we developed to reduce the data 
cubes is described in Kervella & Domiciano de Souza f 20071) . 
After a classical chopping-nodding subtraction of the back- 
ground, the frames were precisely recentered. We then spatially 
resampled each individual frame by a factor four. This oversam- 
pling allows us to take advantage of the atmospheric jitter of 
the images to improve the sampling of the point spread function 
(through a classical dithering technique). To compute the final 
images, we averaged the complete image cubes in all filters, as 
the improvement in angular resolution brought by image selec- 
tion was marginal, while degrading the sensitivity to the faint 
extensions of the nebula. 

2.3. Treatment of image saturation 

The extreme brightness of our two target stars resulted in a sat- 
uration of the VISIR detector on the PSF cores in some of the 
filters. For Aldebaran, some limited saturation occured with the 
PAHl, SIV, and PAH2 filters. For Betelgeuse, all filters in the 
band showed some degree of saturation: low for J7.9 and Nell, 
moderate for SIC, and heavy for PAHl, SIV and PAH2. The 
Ql and Q3 filters did not produce saturation for the two stars. 
The eff'ect of saturation of the VISIR images is two-fold: firstly, 
it causes an underestimation of the total photometry of the satu- 
rated image, and secondly, it creates a number of image artefacts. 

2.3.1. PSF core saturation 

The eff'ect of the photometric bias due to saturation is diff'erent 
for Betelgeuse and Aldebaran. For Betelgeuse, the bias on the 
photometry of the PSF core is not a problem since we are inter- 
ested in the CSE, that, although rather bright, does not saturate 
the detector. For Aldebaran (PSF and photometric calibrator), 
the saturation of the PSF core in the PAHl, SIV and PAH2 fil- 
ters creates a bias on the total flux, and makes the images un- 
usable for the deconvolution of the images of Betelgeuse. To 
overcome this problem, we extrapolated the saturated PSF core 
of Aldebaran using a synthetic PSF profile. It was computed tak- 
ing into account the shape of the telescope pupil (including the 
secondary mirror obscuration), the photometric transmission of 
the atmosphere (with airmass dependence), the transmission of 
the VISIR filters, and the quantum eflftciency of the detector. The 
resulting synthetic PSF images were normalized to the observed 
flux in the images of Aldebaran, using as normalization a ring 
centered on the star and avoiding the saturated area. In practice, 
we integrated the flux over a ring from ^ 0.3'' to ^ 0.9'' from 
the star. The result of the PSF core extrapolation, as well as the 
normalization and integration domains are shown in Fig.|2] The 
normalized synthetic PSF core was used for the extrapolation 
only where its value exceeded that of the original saturated im- 
age. This resulted in a continuous and smooth transition between 
the synthetic extrapolated core and the external part of the PSF. 

2.3.2. Cosmetic artefacts 

Another efl'ect of saturation is the presence of several types of 
cosmetic artefacts, mainly stripes, on the images. The memory 
effect of the detector also causes some ghosts to appear at the 
chopping positions where the star is positioned, but this only re- 
sults in a moderate photometric bias, on the saturated part of 
the PSF, that is unimportant for our analysis. The striping ef- 
fect is more problematic, as it creates a noticeable pattern over 
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Fig. 2. Principle of the restoration of the saturated PSF core 
of the images of Aldebaran. The curves represent radial me- 
dian profiles of the original saturated image (thin red curve), 
the desaturated image (thick green curve) and the synthetic PSF 
used for the central extrapolation (dashed grey curve). The non- 
linearity and saturation limits are represented with dashed and 
dotted grey lines. 




Fig. 3. Stripping of the Betelgeuse average image in the PAHl 
filter caused by heavy saturation (left panel) and corrected ver- 
sion (right). The scale is logarithmic to enhance the faint exten- 
sions of the image. 



the image, that directly affects the photometry of the nebula. 
Fortunately, the removal of this pattern is efficiently obtained 
by subtracting the median value of each line computed over a 
segment near the right and left edges of the frame (Fig.jSj. Some 
saturation-related artefacts are still present close to the star, in 
particular directly above and below, but they remain at an ac- 
ceptable level in all filters. The resulting images of Betelgeuse 
and Aldebaran in all filters are presented in Fig.|4] 

2.4. PSF subtraction and photometric calibration 

The CSE of Betelgeuse is contaminated by the numerous diffrac- 
tion rings from the (unresolved) central star. These rings present 
a significant intensity compared to the nebula, and even domi- 
nate close to the star. It is therefore important to subtract them 
properly. The core of the PSF of Betelgeuse being saturated in 
all N-band filters, the normalization of the PSF images is partic- 
ularly difficult. It is also not possible to use the unsaturated part 
of the images of Betelgeuse (e.g. the external diffraction rings) 
to normalize the PSF image, as the nebula contributes signif- 
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Fig. 4. Average VISIR images of the PSF calibrator Aldebaran (left image of each column) and Betelgeuse (right image of each 
column) in 8 filters from 7.76 to 11.25 fim (left column) and 11.85 to 19.50yum (right column). The color scale is a function of the 
square root of the intensity, scaled to the minimum and maximum values in each image. The field of view is 5.63'' x 5.63'' for all 
images, with North up and East to the left. The images of Aldebaran in the PAHl, SIV and PAH2 filters, as well as all the N band 
images of Betelgeuse include a correction of the PSF core saturation (see Sect.|2.3|for details). 



icantly to the total flux and would cause an overestimation of 
the normalization factor. Due to these constraints, we adopted 
an original approach based on theoretical models of the spectral 
energy distribution of Betelgeuse (the star itself, excluding the 
envelope) and Aldebaran. 

We used the broadband magnitudes of Betelgeuse 
and Aldebaran from the literature to calibrate Castelli & 
Kurucz ( 2003 ) models of their spectral energy distributions (see 
Paper I for details). We used these models to calibrate photo- 
metrically the images of Aldebaran, and derive the ratio of the 



apparent fluxes of the two stars in all our filters (including the 
atmosphere and detector transmissions). We then subtracted the 
scaled images of Aldebaran from the Betelgeuse images. Using 
this procedure, we subtracted the photometric contribution 
from the central star, but we leave the envelope untouched. A 
comparison of an image of Betelgeuse before and after the PSF 
subtraction is presented in Fig.|5] The numerous diff'raction rings 
of Betelgeuse are efficiently subtracted and do not appear on the 
resulting envelope image. The main geometrical characteristics 
of the shell can clearly not be attributed to the PSF subtraction 
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Fig. 5. Image of Betelgeuse in the PAH2 filter before (left) and 
after (right) the subtraction of the PSF image derived from the 
observation of Aldebaran. The grey scale is a function of the 
square root of the intensity between and 0.5 Wm"^yum"^ sr"^ 
The black disk in the center corresponds to the non-linear and 
saturated part of the Betelgeuse image that w^as masked after the 
subtraction. 



as the PSF is circularly symmetric and as the departure from 
symmetry of the shell is visible on the raw images. 

2.5. Image deconvolution 

The VISIR images obtained using the BURST mode are well 
suited for deconvolution, as they are essentially diff'raction lim- 
ited, and mostly insensitive to the seeing. We deconvolved 
the PSF subtracted images of Betelgeuse using the images of 
Aldebaran (corrected for core saturation in the PAHl, SIV and 
PAH2 filters) as the dirty beams, and IRAlQs implementation 
of the Lucy-Richardson (L-R) algorithm (lucy command of the 
stsda^ package). We stopped the L-R deconvolution after a 
number of iterations depending on the wavelength: 5 iterations 
for J7.9 and PAHl, 10 iterations for SIV and PAH2, 20 iterations 
for SIC and Nell, and 40 iterations for the Ql and Q3 filters. 
This variable number of iterations results in an homogeneous 
deconvolved image quality for all filters, while avoiding the cre- 
ation of deconvolution artefacts. This is particularly important at 
the shortest wavelengths (J7.9 and PAHl filters), where the neb- 
ula is faint and the PSF subtraction is delicate. In the Q band, 
there is no saturation, and the PSF subtraction is much easier 
and eflftcient. This allowed us to use a larger number of decon- 
volution steps. The resulting deconvolved images of the nebula 
of Betelgeuse in the six N band and the two Q band filters are 
presented in Fig. [6] 

3. Image analysis 

3.1. Structure of the nebula 

The overall atructure of the nebula appears to be approxi- 
mately triangular, with an inhomogeneous surface brightness. 
Extensions are present up to ^ 3'' from the star to the North, 
South and Northwest (Fig. [7]). An irregular ring-like structure is 
visible at a radius of 0.5 to 1.0'', particularly in the Eastern part 
of the SIV filter image (rectangle in Fig. [7]). For concision, we 
will refer to this structure as "the shell", although we do not have 
evidence from our images of its actual tridimensional shape. As 
shown in Fig.[6j the shell is also detectable in the other N band 

^ http://iraf.noao.edu/ 

^ http://www.stsci.edu/resources/software_hardware/stsdas 
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Fig. 7. Main identified features in the circumstellar nebula of 
Betelgeuse. The background is the image in the SIV filter, with 
an inverted grey scale proportional to the square root of the flux. 



Table 2. Coordinates of the photometric measurement areas 
shown in Fig. [7] in the nebula of Betelgeuse, relatively to the 
star. 
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Notes. The position angle (PA) is counted positively from North to East 
(N=0°, E=90°). The last three lines list the positions of the center of 
the shell measurement window, and its south-eastern and north-western 
corners. 



filters. In the two Q band filters, it is clearly visible in the north- 
western quadrant of the image, at a slightly larger angular radius 
from the star. 

The plumes labeled 1 , 2 and 3 appear to brighten inside the 
shell, and may originate from the star itself. Several bright knots 
are also visible in the nebula (labels A to F, ordered by decreas- 
ing brightness in the SIV filter). Some of these knots may be con- 
nected to the star through the plumes (knots A and C). Additional 
filamentary structures are observed at distances of 2 to y from 
the star (labels 4 and 5), and a large elliptical loop is marginally 
detectable in the SIV, PAH2 and SIC images. The coordinates of 
the photometric measurement windows over the knots and the 
shell are listed in Tabled 
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Fig. 6. PSF subtracted and deconvolved images of the CSE of Betelgeuse for the eight VISIR filters in the N and Q bands. The field 
of view is 5.63'' x 5.63'' for all images, with North up and East to the left. The color scale is a function of the square root of the 
intensity and is normalized to the maximum and minimum intensity in the image. The contour levels correspond to 0.01, 0.05, 0.1, 
0.5 and 1.0 W m"^ jim'^ sr"^ . For the band images, the dark spot close to the center corresponds to the saturated part of the image 
that was ignored in the deconvolution process. 



3.2. Overall spectral energy distribution of the nebula 

We derived the spectral energy distribution of the nebula in our 
eight VISIR filters through an integration of the flux over the 
PSF- subtracted , but non-deconvolved images. We did not use 
the deconvolved images, in order to avoid possible biases due to 
the deconvolution algorithm, and also because we did not need 
the high angular resolution. We integrated the flux of the nebula 
over a circular aperture of 6" in diameter. The central part of the 
nebula is lost during the PSF subtraction, due to the saturation 
of the core of the images of Betelgeuse (except in the Q band). 
To avoid a bias on the integrated photometry, we extrapolated 
the nebular flux in the central masked disk (0.4" in radius) us- 
ing a Gaussian fit over the rest of the PSF- subtracted image of 
the envelope. This resulted in practice in a marginal additional 
contribution for all filters (< 3%), except in the J 8. 9 and PAHl 
filters 20% and 10%, respectively). 

The estimation of the uncertainty in the determination of the 
absolute flux of Betelgeuse is difficult, due to the treatment of 
saturation that we had to apply to the PSF reference images. 
Thanks to the high brightness of the two stars, the statistical dis- 
persion of the photometry is negligible, compared to the system- 
atics. We estimate that the systematic uncertainty on our mea- 
surement of the irradiance of Aldebaran is ^ 20%, and that it 
directly applies to the photometry of the nebula of Betelgeuse. 
In Fig.[8j we show the total irradiance of the circumstellar neb- 
ula as a function of wavelength. The systematic uncertainty on 
our model of the irradiance of Betelgeuse itself (the star alone, 
for the PSF subtraction) is estimated also to 20%, due essen- 
tially to the variability of the star. Overall, this results in ^ 30% 
error bars for the N band measurements (Betelgeuse model un- 
certainty plus saturation correction uncertainty), and ^ 20% in 
the Q band (Betelgeuse model uncertainty alone). 



Fig.[8]also shows for comparison the spectrum of Betelgeuse 
obtained using ISO's (Kessler et al. 1996) SWS instrument (de 
Graauw et al. |1996K retrieved from the ISO Data Archiv^ The 
overall shape of the spectrum is comparable to our VISIR pho- 
tometric data points, although slightly brighter. This systematic 
diff'erence at a 1 - 2cr level could be explained by the relatively 
large aperture diameter of the SWS spectrograph (14 to 27", vs 
6" for our flux integration disk), or by the temporal variability of 
Betelgeuse (the ISO spectrum was obtained on 1997-10-08). The 
absolute photometric uncertainty of SWS low-resolution spectra 
is approximately 5-10%, and may also explain part of the ob- 
served difl'erence. 



3.3. Spectral energy distribution of individual features 

We measured the surface brightness of the difl'erent identified 
features of the nebula on the PSF- subtracted and deconvolved 
images of Betelgeuse. We prefered the deconvolved ima ges for 
this task for two reasons. Firstly, McNeil & Moody (I2005t 
showed that the Lucy-Richardson algorithm preserves the pho- 
tometric accuracy relatively well compared to other classical 
deconvolution algorithms. Secondly, our images cover a very 
broad range of wavelengths, with more than a factor 2.5 between 
our extreme wavelengths (J7.9 and Q3 filters). The diffraction- 
limited angular resolution therefore also varies by this factor, 
and the surface photometry over a given aperture size would re- 
sult in very diff'erent actual measurement areas on the sky de- 
pending on the fflter (the adjacent regions to the aperture would 
contribute diff'erently depending on the angular resolution). For 
these reasons, we prefered to use the deconvolved images for the 
photometry presented in this Section. The azimuthally-averaged 



http://iso.esac.esa.int/ida/, observation ID 69201980 
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Fig. 8. Integrated irradiance of Betelgeuse plus the nebula 
(black dots), superimposed on the model spectrum of the star 
alone (dashed curve), compared to the ISO SWS spectrum of 
Betelgeuse. The horizontal segments attached to each point in- 
dicate the spectral width at half transmission of the VISIR filters. 
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Fig. 9. Azimuthally-averaged surface brightness profile of the 
CSE of Betelgeuse, computed from the PSF-subtracted and de- 
convolved images. 



surface brightness of the CSE of Betelgeuse is represented in 
Fig.|9]as a function of the angular separation from the star. The 
approximate radius of the shell is shown for convenience. 

3.3.1. The shell 

The approximate location of this ring-like structure, that we ten- 
tatively interpret as a shell, is represented with a dashed circle 
in Fig. |7] Although its overall shape is circular, it is relatively 
irregular in its western section, whereas its eastern part shows a 
mostly continuous arc in the N band. In the Q band, the western 
section is more clearly defined (Fig. [6]). We nevertheless chose 
the eastern section of the shell to estimate its surface brightness 
(rectangular aperture drawn in Fig. |7]), as it is visible in all N 
band filters. The coordinates of the comers and center of the 
measurement window are listed in Table |2] The resulting pho- 
tometry is presented in Fig. fTOl The shell appears particularly 
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Fig. 10. Average surface brightness over a rectangular area lo- 
cated on the Eastern section of the shell (defined in Fig. [7]). The 
vertical error bars correspond to the standard deviation of the 
surface brightness over the measurement area, and the horizon- 
tal error bars indicate the spectral width at half transmission of 
the VISIR filters. 



bright in the SIV, PAH2 and SiC filters which is not systemati- 
cally the case closer to the star. Another interesting property is 
that it appears brighter in the Q3 band (10.7 ± 1 . 1 mW/m^/^m/sr) 
than in the Ql band (3.2 ± 1.2 mW/m^A/m/sr). This spectro- 
photometric pattern is compatible with the presence of typical 
0-rich dust, such as silicate-type materials. In the average ra- 
dial profile of the CSE shown in Fig. [9] the photometric contri- 
bution of the shell is visible in the SIV, PAH2 and SIC filters 
as a "bump" in the azimuthaly-averaged profile. It is also eas- 
ily observable in the Q band filters, thanks to its higher contrast 
compared to the nebular background. 

3.3.2. Nebular knots 

The average surface brightness of the six knots identified in the 
nebula are presented in Fig.fTT] We measured the surface bright- 
ness over a circular aperture of 0.38'' in diameter, a size that cor- 
responds to the angular resolution of VISIR around 13 yum. This 
wavelength is around the center of our full wavelength range, 
and combined with the image deconvolution we applied, this rel- 
atively large aperture results in surface brightness measurements 
that are largely insensitive to the changing angular resolution 
with wavelength. As discussed in Sect. |3.2| we estimate the sys- 
tematic uncertainty on the photometry of the nebula to 20%. The 
error bars represented in Fig. [TT] correspond to the dispersion 
of the flux on each of the measurement apertures, quadratically 
added to the systematic uncertainty. 



4. Discussion 

4A. Evolutionary time scales 

An important aspect to consider in any study of Betelgeuse is 
its variability with time. This makes the comparison of difl'er- 
ent studies obtained at difl'erent epochs more difficult, but it is 
also a powerful means to sample the physics of the star, espe- 
cially since the advent of realistic 3D hydrodynamical simula- 
tion of supergiant stars (Frevtag 120021 Freytag & Hofner [20081 
Chiavassa et al. .2010i) . Betelgeuse exhibits a variability over two 
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Fig. 11. Average surface brightness of Features A, B, C (top) 
and D, E, F (bottom) as defined in Fig. [7] measured over a disk 
with a diameter of 0.38''. The error bars correspond to the stan- 
dard deviation of the surface brightness over the measurement 
area, added quadratically to the 20% systematic caHbration un- 
certainty (for the N band only, see text for details). 



periods of approximately 200 and 2100 days (Stothers "2010), 
by a few tenths of magnitude in visible photometry, and a few 
kms"^ in radial velocity. From spectroscopic observations in 
the visible domain. Gray ( 2008 ) report a characteristic evolution 
timescale of the velocity variations at the surface of 400 days, 
and a range of velocities of 9 km s"^ 

Convective-related surface structures are characterized by 
two characteristic time scales (Chiavassa et al. 1201 11 : (1) large 
convective cells, with a size comparable to the stellar radius, that 
evolve on time scale of years; these convective cells are visible in 
the infrared. On the top of these cells, there are small convective 
structures evolving on time scale of months. (2) In the optical re- 
gion, large convective cells are still visible but they are hidden by 
numerous short-lived (a few weeks to a few months) small-scale 
structures (about 5 to 10% of the stellar radius). An observa- 
tional evidence of this long timescale is the probable persistence 
of the southwestern plume reported in Paper I, that probably cor- 
responds to the bright spot noticed more than 10 years before 
in HST images of the chromosphere of Betelgeuse (Uitenbroek 
et al.|W98, see also Gilliland & Dupree 1996 ). 

At higher altitudes in the CSE, we can derive a rough esti- 
mation of the evolutionary time scale from the ratio of the lin- 



ear distances involved divided by the displacement velocity of 
the gas. The velocity of the CO gas was estimated recently by 
Smith et al. (2009 ) from spatially resolved spectroscopy 8it A = 
4.6yum. They observed typical velocities of 12kms"\ in very 
good agreement with the slow gas component detected by Bemat 
et al. (119791) for a temperature of ^ 200 K, who additionally re- 
port a faster expanding batch of gas (-hi 8 km s"^) at T^^ 70 K. To 
convert angular velocities in linear velocities, we use the revised 
Hipparcos parallax of Betelgeuse from van Leeuwen (12007b , 
71 = 6.55 ± 0.83 mas, a value compatible with the recent results 
obtained by Harper et al. (2008 ). Considering this distance of 
152 ± 20 pc, the photosphere (~44 mas in the K band, Perrin 
et al. 2004) has a linear radius of 3.4 AU, and the shell (angular 
radius of 0.9'') is located at ^ 134 AU. Considering a gas veloc- 
ity of 10-20 km s"\ the gas will take ^ 30-60 years to migrate 
from the photosphere to the inner radius of the dust shell. 

In addition to a physical displacement of matter, the change 
in illumination of the nebula by the evolving photosphere of the 
star can result in an apparent morphological change of the neb- 
ula, as proposed by Skinner et al. (119971 ). This evolution will 
happen with a typical timescale comparable to the diff'erent con- 
vective structures present on the star's surface (months to years). 
Using the available instrumentation, it therefore appears possi- 
ble to monitor the evolution of the photosphere of Betelgeuse 
over a few months using interferometry, and of its CSE over a 
few years using AO for the inner gaseous envelope, and thermal 
infrared imaging for the dusty envelope. 



4.2. Comparison with NACO images 

A color composite image of the present VISIR observations 
and the near-infrared observations in the JHK bands (Paper I, 

The two 



shown in the central disk) is presented in Fig. 12 
imaging data sets are represented to scale, the central white 
spot corresponding to the photosphere of the star. One of the 
key questions we wish to address in our observing program of 
Betelgeuse is whether CSE asymmetries (plumes in the near- 
IR, nebular features in the thermal IR) are linked to the convec- 
tive cells imaged by interferometry at near-infrared wavelengths 
(Haubois et al. [2o69l Ohnaka et al. 2009, Ohnaka et al. [2OTT1 
see also Chiavassa et al. 120091 ). From the arguments presented in 
Sect. 4.1 it is justified to directly compare the observations ob- 
tained with NACO (obtained in January 2009) with the present 
VISIR images (obtained in November 2010), as the evolutionary 
time scale of the dusty envelope is significantly longer than the 
two years separating the two data sets. 

Considering their direction and opening angle, the southern 
and northwestern near-infrared plumes extending from the pho- 
tosphere correspond well with the thermal infrared filamentary 
structures labeled #2 and 3 in Fig.|7] This is an indication that the 
molecular flow from the photosphere of the star extends through 
and feeds the dust shell. The faint loop observed at larger dis- 
tance from the star (#6 in Fig. [7]) could correspond to a past 
episode of enhanced mass loss. From the spectral energy dis- 
tribution of the nebula, it is clear that the shell contains a sig- 
nificant ammount of 0-rich dust (silicates, alumina,...), which 
is expected for Betelgeuse (see e.g. Verhoelst et al. .2009 j . This 
type of dust has strong signatures around 9-lOyum, that are vis- 
ible in the ISO spectrum presented in Fig. [8] as well as in our 
photometry of the nebula (Fig. [8]). It is expected to form at a few 
0.1" above the stellar surface, which also corresponds well to 
the radius of the shell measured in the VISIR images. 
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4.3. Previous observations of Betelgeuse in thermal and far 
infrared 

The dusty environment of Betelgeuse has been studied at thermal 
infrared wavelengths using several complementary techniques: 
interferometry (at a few tens of milliarc second scale), large scale 
imaging (a few arcsecond scale) and spectroscopy. 

Interferometric observations at llyum with the Infrared 
Spatial Interferometer (ISI) allowed Bester et al. (119911 to mea- 
sure the inner radius of the dust shell at 0.90 ± 0.05", which 
was confirmed by Danchi et al. J1994b . This radius is in very 
good agreement with the shell we observe with VISIR (Fig.[7|. 
Although circular symmetry is always hypothesized, it appears 
in our images that this shell is not complete nor uniform. Townes 
et al. (2009) reported a decrease of the angular diameter of 
Betelgeuse by 15% between 1993 and 2009, from interfero- 
metric measurements (also with ISI), using two-telescope (up 
to 2002) and 3-telescope (2002-2009) configurations. While the 
changing fringe visibility may indeed be related to a change 
in the photospheric size of the star, another explanation could 
be that the complex structure of the CSE aff'ected the measure- 
ment. Its high brightness could induce a significant variation of 
the interferometric visibility with baseline length and orienta- 
tion, that would result in a deviation from the single star model, 
and consequently a bias on the derived angular diameters. Tatebe 
et al. (120071) identified an elongation of the disk of Betelgeuse 
at 11.15yum using ISI with 3 telescopes (their Fig. 2), approxi- 
mately in the North-South direction. To the North of Betelgeuse, 
our VISIR images show a bright, jet-like structure (#1 in Fig. [7]). 
It is clearly visible in the 11.25yum VISIR image (Fig.|6]), and 



^^^^^^^^^^H Fig. 12. Color composite image of the CSE of 

^^^^^^^^^^H Betelgeuse. The NACO observations obtained 

^^^^^^^^^^H in January 2009 (Paper I) are reproduced in the 

^^^^^^^^^^H central disk, with J coded as blue, H as green 

^^^^^^^^^^H and K as red. The apparent angular extension of 

^^^^^^^^^^H the near-infrared photosphere (^ld ^ 44 mas) 

^^^^^^^^^^H is represented with a red circle at the center of 

^^^^^^^^^^H the image. The field of view is 5.63'' x 5.63'', 

^^^^^^^^^^H with North up and East to the left. The thermal 

^^^^^^^^^^H infrared VISIR images reported in the present 

^^^^^^^^^^^H work are presented with PAHl coded as yellow, 

^^^^^^^^^^H SIV as orange, and Ql as red. 

its location well within the field of view of ISI may explain the 
measured apparent elongation. 

An original intermediate approach between long base- 
line interferometry and direct imaging was chosen by Hinz 
et al. ( 1998 ). Using two mirrors of the Multiple Mirror Telescope 
separated by a baseline of 5 m and nulling interferometry, these 
authors observed the CSE of Betelgeuse at lOyum after cancel- 
lation of the flux from the star itself. They determined a trian- 
gular shape for the nebula that we also observe in o ur ima ges, 
with the same orientation and extension. Harper et al. ("2009) ob- 
served Betelgeuse using the TEXES spectrograph installed on 
NASA's IRTF, and identified narrow iron emission lines at 17.94 
and 24.52 yum, that form close to the stellar surface (within a few 
stellar radii). Such emission could explain at least part of the un- 
resolved emission we observe in the Ql and Q3 VISIR images 
at the star's location (Fig. [6]) 

Few direct imaging observations of Betelgeuse in the ther- 
mal infrared domain are present in the literature. Skinner et 
al. (1997 ) obtained images of Betelgeuse at 8.2, 9.7 and 12.5 jim 
using UKIRT. They resolved the circumstellar nebula, and no- 
ticed a significant evolution over a period of one year. They 
attribute this evolution to a change in illumination of the neb- 
ula, caused by the presence of variable spots on the star's sur- 
face. From a model that includes also radio observations with 
the VLA (see also the preceding modeling work by Skinner & 
Whitmore 1987 ), they determine an internal radius of the dust 
shell of 0.5'', which is slightly smaller than the shell we observe 
with VISIR. Direct imaging of Betelgeuse in the thermal infrared 
was also obtained by Rinehart et al. (1998 ) at 11.7 and 17.9yum 
using the Hale 5 m telescope. From a model fit to their images, 
they derive an inner dust shell radius of 1.0 ± O.T', also in excel- 
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lent agreement with our observations. A very similar radius for 
dust formation (^^33 = 0.73'') was also proposed by Harper et 
al. ( 2001 ) from a model of the extended atmosphere of the star. 

At a much larger spatial scale, the interface between the wind 
of Betelgeuse and the interstellar medium was recently observed 
in the far infrared (65 to 160 yum) by Ueta et al. (120081) using the 
AKARI satellite. This bow shock is located at a radius of 4.8' 
from the star, corresponding to ^ 100 times the extension of the 
dusty envelope we observe with VISIR. 

4.4. A stellar companion ? 



Some interesting structures were detected by Smith et al. (120091) 
^ 1" west of the star, which corresponds to the position 
of the bright knot A (Fig. [7|. Its spectral energy distribution 
(Fig. [TT]) indicates that it contains 0-rich dust. Its compact as- 
pect and relatively high brightness are however intriguing, and 
it is not excluded that it may correspond to a stellar compan- 
ion of Betelgeuse. Such a companion would be in a position 
to accrete material lost by the supergiant, and therefore present 
dust and molecular emission features. Considering the mass of 
Betelgeuse 10 - 20 M©), its a priori probability to be part 
of a binary or multiple system is relatively high 70%, see 
e.g. Kouwenhoven et al. 12007b . It should be noted that knot A 
is located outside the field of view of the NACO observations 
reported in Paper I, and we cannot confirm this hypothesis with 
these observations. 

The presence of two companio ns in o rbit around Betelgeuse 
was proposed by Karovska et al. (11986b based on speckle in- 
terferometric observations obtained at epoch 1983.88, but they 
were not confirmed. One of the sources detected by these au- 
thors was found at a position angle of 278 ± 5° and a radius of 
0.51 ± 0.01" from Betelgeuse, i.e. closer but at a similar position 
angle compared to our knot A (Table[2]and Fig. [7]). If we consider 
that it is the same source observed at two epochs, its angular ve- 
locity projected on the plane of the sky would be 33masyr"^ 
(VISIR epoch is 2010.86). At the distance of Betelgeuse, this 
angular velocity corresponds to a projected linear velocity of 
^ 24kms"^ relatively to Betelgeuse. This figure is compara- 
ble to the typical velocities derived from CO lines, considering 
the diff'erent projection eff'ect (spectroscopy measures the veloc- 
ity perpendicular to the plane of the sky, while imaging gives 
measures it in the plane of the sky). However, this apparent dis- 
placement appears too slow for the orbital motion of a close-in 
companion of Betelgeuse, that would have a typical period of a 
few years at a separation of 1". The possibility nevertheless re- 
mains that knot A corresponds to a relatively distant companion 
(thus with a slower orbital motion), whose apparent proximity to 
the star would be caused by projection eff'ects. 

While this hypothesis cannot be formally excluded from our 
data, the possibility that knot A is associated with the presence 
of a stellar companion of Betelgeuse appears unlikely. The most 
probable hypothesis is that it is simply an overdensity of the dust 
in the CSE, that ma y have been observed in visible scattered light 
by Karovska et al. (11986b . 

5. Conclusion 

Our VISIR observations of the CSE of Betelgeuse show the pres- 
ence of a complex and bright nebula. The observed brightening 
of the nebula for wavelengths greater than ^^9-10 yum points at the 
presence of 0-rich dust, such as silicates or alumina. The shell 
that we observe at a radius of 0.9" presents a similar spectral en- 
ergy distribution as the rest of the nebula, and could correspond 



to the condensation radius of certain dust species, or to a past 
episode of enhanced mass loss. A promising correspondance is 
observed between the plumes reported in Paper I from NACO 
observations and the VISIR nebular features. The knots and fil- 
amentary extensions of the nebula observed at larger distances 
from Betelgeuse appear to correspond to inhomogeneities in the 
mass lost by the star in the recent past, probably within the last 
few centuries. Further observations are expected to clarify the 
nature and composition of the nebular features identified in our 
images, using spatially resolved spectroscopy of the CSE. 
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